Geometric phases induced in auxiliary qubits by many-body systems near its critical 

points 
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The geometric phase induced in an auxiliary qubit by a many-body system is calculated and 
discussed. Two kinds of coupling between the auxiliary qubit and the many-body system are 
considered, which lead to dephasing and dissipation in the qubit, respectively. As an example, 
we consider the XY spin-chain dephasingly couple to a qubit, the geometric phase induced in the 
qubit is presented and discussed. The results show that the geometric phase might be used to 
signal the critical points of the many-body system, and it tends to zero with the parameters of the 
many-body system going away from the critical points. 
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I. INTRODUCTION 



Motivated by the fact that all realistic quantum sys- 
tems are coupled to their surrounding environments, 
many researchers have paid their attention on the study 
of geometric phase in open systems since 1 980s [H 0,13, Q, 
0, IE 0, 0, 0, 51 [Q, dlH 51 • From the perspective 
of application, the use of geometric phase in the imple- 
mentation of fault-tolerant quantum gates [H, [l?], [H, HI] 
also requires the study of geometric phase in open sys- 
tems, because the system which needs to be manipulated 
may decoher from a quantum superposition into statis- 
tical mixtures due to coupling to its environment. De- 
spite of extensive studies that have been performed in 
this direction, investigation into the effect of correlation 
among the environmental systems on the induced geo- 
metric phase in open systems is lacking. 

Geometric phases have been shown to associate with a 
variety of condensed matter phenomena [If], Hl|, HU, [H| . 
Nevertheless, their connection to quantum phase transi- 
tions has only been given recently in [H [25], [26| , where 
it has been shown that the Berry's phase can be used to 
signal the critical points of the spin chain 24] , the criti- 
cal exponents were evaluated from the scaling behavior 
of geometric phases 25], and the geometric phase can be 
considered as a topological test to reveal quantum phase 
transitions [26]. In these works, the geometric phase an- 
alyzed is acquired by the ground state (or the low lying 
excited states) of the many-body system. This gives rise 
to the following question, whether the geometric phase 
induced in an auxiliary system by the many-body sys- 
tem has connection to the critical points of the many- 
body system. In this paper, we shall try to answer this 
question by examining the geometric phase induced in a 
qubit by the many-body system. This makes our study 
different from that in [24], [25], [26|. Two kinds of cou- 
pling between the auxiliary qubit and the many-body 
system are considered. The first kind of coupling com- 
mutes with the free Hamiltonian of the auxiliary qubit, 
leading to dephasing in the qubit. Whereas the second 
one commutes with the free Hamiltonian of the many- 



body system, which causes energy loss (dissipation) of 
the qubit. 

This paper is organized as follows. In Sec. II we present 
a general formulism to calculate and analyze the reduced 
density matrix of the auxiliary qubit, and establish a 
connection of the reduced density matrix to the criti- 
cal points. As an example, we calculate the geometric 
phase induced in the auxiliary qubit taking the XY spin 
chain as the many-body system in Sec. III. And finally we 
conclude our results in Sec. IV. 



II. GENERAL 

In this section, we will give a general formalism to an- 
alyze the reduced density matrix of a qubit coupled to 
a quantum many-body system, and establish a connec- 
tion between the reduced density matrix of the qubit and 
the criticality in the quantum system. We restrict our- 
self to consider the following qubit to many-body system 
couplings. First we consider the case where the coupling 
conserves the energy of the qubit. Then we analyze the 
case where the energy of the qubit does not conserve, 
but the energy of the quantum many-body system con- 
serves. The first case corresponds to dephasing in the 
qubit, while the second kind of coupling results in dissi- 
pation in the qubit. 

Consider a qubit coupled to a quantum many-body 
system. The Hamiltonian that governs the evolution of 
the whole system may have the form 



H = Hi 



-ff m (A) + H, 



(1) 



where Hi = /j| f)("f |, describes the free Hamiltonian of 
the qubit, H m (X) stands for the free Hamiltonian of the 
quantum many-body system, and Hi = gHu ®\ |) (| | 
represents the coupling between them. Hm is an arbi- 
trary operator of the many-body system. It is clear that 
[Hi, Hi] — 0, therefore the energy of the qubit conserves. 
The quantum system described by H m (X) undergoes a 
quantum phase transition for parameter A = A c . It is 
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easy to show that the time evolution operator for the 
whole system may be written as 

C/(*) = C/T(*)lt)(tl + ^(t)|l)ai, (2) 
with E/f (f) and Ui(t) satisfying 

ih^U^it) = H^ tl U u (t), (3) 

where 

H yi =H m + {t ±f) + (f ±\)H M - (4) 

|) and | I) are eigenstates of a z for the qubit. Having 
these expressions, we now show that the off-diagonal ele- 
ments of the reduced density matrix of the qubit change 
drastically in the vicinity of critical points. To this end, 
we assume that the qubit and the quantum many-body 
system are initially independent, such that we take the 
following product state as the initial state of the whole 
system 

|V(0)> = |^(0)) ® |G m >, (5) 

where \G m )) represent the ground state of H m and 
iV'i(O)) — cos#| j) + sin#| I). By a standard calcula- 
tion, one obtains the reduced density matrix of the qubit 

_ / cos 2 6» cos9sm9F(t)\ , . 

P h ~ \cos8sm9F*{t) sin 2 9 )' 

Here, F(t) is defined by 

F(t) = (G m \e-i H ^\G m ), (7) 

this is the survival probability of the ground state of the 
quantum system under the action of the Hamiltonian . 
The leading term F(t) ~ (G m |Gj„) of this equation rep- 
resents the overlap function between two ground states 
\G m ) and |G]„), where IGj^) denotes the ground state of 
Hp This overlapping function was shown [27| to take 
extremal values in the vicinity of critical points. In fact, 
F(t) represents the Loschmidt echo which goes exponen- 
tially to zero as the many-body system approaches the 
regions of criticality [28| . This give us the behavior of F(t) 
in the vicinity of critical points A c as F(t) ~ |A — X c \ v , 
leading to a dramatic change in the induced geometric 
phase near this point A c . We would like to notice that 
the above discussions hold for a general coupling between 
the spin and the quantum system, provided the coupling 
is weak. 

Next we turn to study the case when the energy of the 
qubit does not conserve, but the energy of the many-body 
system conserves instead. This implies that [Hi, Hi] ^ 0, 
but [Hi, H m ] — 0. Without loss of generality, we consider 
the following Hamiltonian 

H = Hi+H m (X)+Hi, 
Hi = Aa x , Hi = (g z cr z + g y <J y ) ® H m (X). (8) 



<Tj(« = x,y,z) are Pauli operators of the qubit. The in- 
teraction Hamiltonian Hi can be rewritten as 

H i = ge i ^a z e- i ^, (9) 

with cos(27) = g z /g, and g — y/ g z + g 2 z . The fact that 
the interaction Hamiltonian Hi commutes with H m en- 
ables us to write the time evolution operator of the whole 
system as 

N 

U(t)=Y,Un(t)\EnW)(EnWl i 1 ®) 

71=1 

where {|-E„(A))} are eigenstates of H m with correspond- 
ing eigenenergies {E n = E n (X)}, and we assume that 
they are nondegenerate. It is easy to show that 

9 9 
U n (t) — e~ nt {cos — cos Q n t — io x cos — sin £l n t 



+ i sin — cos Q, n t[(cos 2j + sin 2j)a z 

+ (sin 27 - cos27)cr y ]}. (11) 

Here, cos(9„ = , and ft n = Jg 2 E^+A 2 . 

Quantum phase transition theory tells us that the ground 
state energy Eq(X) behaves drastically in the vicinity of 
the critical point A c , this would reflect in Uo(t) that is 
a function of Eq(X). Having established this linkage, we 
claim that the final state of the qubit 

P i (t)=^|c n | 2 [/ n (t)pi(0)C/t Wj (12) 

n 

and the geometric phase acquired by this state can sig- 
nal quantum critical points in the many-body system, 
provided the initial state of the many-body system is 
a superposition of {\E n )}, i.e., |^ m (0)) = Y,n=o c n\ E n) 
(co 7^ 0), and at least two c„ including Co are not zero. We 
would like to emphasize that the geometric phase in this 
case can not signal the critical points, if the many-body 
system is initially in the ground state with probability 1 . 
This is because the qubit-system interaction could not 
excite the many-body system in this situation. In other 
words, if the couplings between the qubit and the many- 
body system commute with the free Hamiltonian of the 
many-body system, the many-body system would remain 
in its initial state if the initial state is one of the eigenval- 
ues of the free Hamiltonian H rn . Thus, the many-body 
system would make no effect on the qubit during the dy- 
namics, and consequently the geometric phase acquired 
by the qubit could not signal the critical points of the 
many-body system. 

III. GEOMETRIC PHASES INDUCED BY THE 
XY SPIN CHAIN 

In this section, we will present an example to illus- 
trate the claim. Taken a spin-chain described by the one- 



3 



dimensional XY model as the quantum many-body sys- 
tem, the Hamiltonian for the whole system (spin+chain) 
may be given by 

H = Hi + H m + H i: (13) 

where 
Hi = ps z , 

JV 

H m = -2^((l+ 7 ) s f s f +1 + (l- 7 ) s f s f +1 + A S f), 

Hi = ig^s'sf. (14) 
i=i 

Here s denotes spin operator of the qubit which cou- 
ples to the chain spins s; (I — l,...,N) located at 
the lattice site I. The spins in the chain are coupled 
to the qubit through a constant g. The time evolu- 
tion operator U(t) for the whole system may be writ- 
ten as U(t) — Uj{t)\j){j\. It i s eas Y to show 
that Uj(t)(j =T,J.) satisfy ih^Uj(t) = H 3 Uj(t) with 

Hj = -E^iC^offff+i + 1 T La i°i+x + A j°f); where 
Xj = A ± g. If the auxiliary qubits are initially in state 
\j), the dynamics and statistical properties of the spin 
chain would be govern by Hj, it takes the same form as 
H m but with modified field strengths Xj. The Hamilto- 
nian Hj can be diagonalized by a standard procedure [29j 

to be Hj = T,k u iA'0j,kV3,k ~ |), where VjAv}^) are 
the annihilation (creation) operators of the fermionic 

modes with frequency uij.k — J e 2 k + j 2 sin 2 Sj t k = 



(cos 2# - Aj), k = -N/2, -N/2 + 1, ...,N/2 - 1. The 
fermionic operator rjj k was defined by the Bogoliubov 

transformation as, rjj^k = dfc cos — id)_ k sin -4^, where 
dk = —j= E; a / ex P( — i2irlk/N), and the mixing angle 0^ 
was defined by cosBj^ — Sj,k/<^j,k- Fermionic operators 
ai were connected with the spin operators by the Jordan- 
Wigner transformation via ai — (Y\ m<l <T m)( <7 i + 
The operators rjj^ parameterized by j clearly do not 
commute with each other, this will leads to nonzero ge- 
ometric phase in the auxiliary qubits as shown later on. 
Before going on to calculate the reduced density ma- 
trix, we present a discussion on the diagonalization of 
Hj. For a chain with periodic boundary condition, i.e., 
o\ = crjV) boundary terms Hi, oun ~ [(ajyOi + 7ajyai) + 
h.c.][exp(iirM) + 1] have to be taken into acount [29l |3C| . 
In this paper, we would work with Hf, un ignored, because 
we are interested in finding a link between the critical- 
ity of the chain and the entanglement in the auxiliary 
qubits. In fact, the boundary effect can be ignored pro- 
vided N — > oo, this is exactly the situation we consider 
in this paper. 

Having given an initial product (separable) state of the 
total system, |^(0)) = \ipi(0)) ®|<^ m (0)), we can obtain 
the reduced density matrix for the auxiliary qubits as 
pi{t) = Tr m [U{t)\i>(0)){ip(0)\U^ (t)], it may be formally 
written in the form 

id 

A straightforward but somewhat tedious calculation 
shows that pij(t) = Pij(t — O)-Fy(t), with 



I 

Fij{t) = JJe^"'.*-^.^^! - (1 - e^. fct )sin 2 ° k - ° hk - (1 - e -*"***) sin 2 ° k - 6]M 

k 

+ (1 - e **.»«)(l - e -**.»«)[sin ° k ~ 2 ° l ' k sin ° k cos ° l ' k ~ ^ k }}, (16) 



where cos 6i. 



cos(2irk/N)-\ 



To derive 



• s /(cos(27rfc/JV)^A) 2 +7 2 sin 2 (2wk/N) ' 

this result, the spin chain was assumed to be initially in 
the ground state of H m . 

With these expressions, we now turn to study the geo- 
metric phase of the open system. For an open system, the 
state in general is not pure and the evolution of the sys- 



tem is not unitary. For non-unitary evolution as shown 
in Eq. lfTB]) . the geometric phase can be calculated as fol- 
lows. First, solve the eigenvalue problem for the reduced 
density matrix pi(t) and obtain its eigenvalues Sk(t) as 
well as the corresponding eigenvectors \ipk(t)}; Secondly, 
substitute and \ipk(t)) into 



$ 9 = arg(]T v / £fc (0)e fc (T)(^(0)|V^(T)) e ^o T ^( t )l a / a *l^( t )> dt ). (17) 

k 

I 



Here, $ 9 is the geometric phase for the system undergo- ing non-unitary evolution [3l|, T is the total evolution 
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FIG. 1: The geometric phase as a function of A. The figure 
was plotted for TV = 1245 sites, g = 0.1 and (a)7 = 0.2,(b)7 = 
0.5,(c) 7 = 0.8,(d) 7 = 1. 



time. The geometric phase Eq. p7|) is gauge invari- 
ant and can be reduced to the well-known results in the 
unitary evolution. It is experimentally testable. The 
geometric phase factor defined by Eq. lfTT)) may be un- 
derstood as a weighted sum over the phase factors per- 
taining to the eigenstates of the reduced density matrix, 
thus the detail of analytical expression for the geometric 
phase would depend on the digitalization of the reduced 
density matrix Eg . (fir?)) . The eigenvalues of the reduced 
density matrix are readily calculated 



e±(t) = ^±^l-4piiP22+4| Pl2 | 2 . (18) 

In order to calculate the geometric phase, we only need 
the eigenvalue £+(t) and its corresponding eigenvector 



this model in the limits of |A| — > oo, |A| = 1 and A = 0. 
The ground state of the spin-chain approaches a prod- 
uct of spins pointing the positive/negative z direction in 
the |A| — > oo limit, whereas the ground state in the limit 
A = is doubly degenerate under the global spin flip by 



n 



N 

1=1 u l 



At |A| 



1, a fundamental transition in the 



ground state occurs, the symmetry under the global spin 
flip breaks at this point and the chain develops a nonzero 
magnetization (a x ) ^ which increases with A growing. 
The above mentioned properties of the ground state are 
reflected in the geometric phase as shown in figure [l]-(d). 
In the limit |A| — * oo, 9j t k = 6k = 7r), this results 
in Q g = 0. In fact when |A| > 4, $ 9 approaches zero 
very well. Around |A| — > 1, the geometric phase changes 
drastically, this can be interpreted as the sensitivity of 
the spin-chain ground state to perturbations from the 
qubit-chain coupling at these points. The difference be- 
tween 7 = and 7 = 1 is that Fij(t) = 1 for 7 = 0, 
but it does not hold for 7=1. This is the reason why 
the geometric phase takes different values at these criti- 
cal points, as figure [1] shows. The ground state of the XY 
model is really complicated with many different regime 
of behavior [33, these are reflected in sharp changes in 
the geometric phase across the line |A| = 1 regardless of 
7 (as shown in figure l-(a),(b) and (c)), indicating the 
change in the ground state of the spin-chain from para- 
magnetic phase to the others. Instead of discussing the 
scaling behavior of the induced geometric phase, we con- 
sider the scaling property of F^j (t) given in Eq. (|16p in 
the vicinity of the critical points. Noticing that i 7 ^ (t) is 
a function of (9^ — 6j,k) (j =T) l)i the scaling behavior of 
Fij(t) may be characterized by S^(X, 7) = Efcli(^r) 2 > 



and S^-(A,7) 

S&(|A| = 1, 7 ) 
in the vicinity of criticality for the one-dimensional XY 
Model. 



It was shown that 
N 2 h 2 and S$f(\\\ < 1,7 = 0) ~ N 2 



since 



MO) 



cos — e 
2 



T) 



•sin 77 1 I), 



(19) 



0. 



Here 
and 



i£li ^/V / |pi2| 2 + i|pii-p 22 p, and tan. 
iTO(pi2)/Re(/?i2), i.e., the ratio of the imaginary 
part to the real part of p\2- Having these results, we 
obtain the geometric phase, 



rg(VM0)e+(T)( £+ (0)| £+ (T)) f 



Jo 



cos %<pdt 



)• 

(20) 

We would like to notice that the geometric phase vanishes 
with 7 = because Fij(t) = 1 at this point. Fij(t) = 
1 results in <f> = in Eq.([20|). this together with that 
(£+(0)|e+(T)) is real yield zero geometric phase at this 
critical point. 

Now we turn to study the criticality in the transverse 
Ising model(7 = 1 in the XY model). The ground state 
structure of this model change drastically as the param- 
eter 7 is varied. We first summarize the ground states of 



IV. CONCLUSION 

In conclusion, by discussing the geometric phase in 
the auxiliary qubit coupled to the many-body systems, 
the relation between the geometric phase induced in the 
qubit and the critical points of the many-body system 
was established. The induced geometric phase in the 
qubit change drastically at the critical points, this is due 
to the sensitivity of the many-body systum to param- 
eter changes near its critical points. The relation not 
only provides an efficient theoretical tool to study quan- 
tum phase transitions, but also proposes a method to 
measure the critical points in experiments 33]. The lim- 
itation of this discussion is that the coupling between 
the qubit and the quantum system is assumed weak, and 
as we have shown, the geometric phase could not reflect 
the critical points of the quantum system that is initially 
in its ground state with qubit-system coupling satisfying 
[H m , Hi] = 0. 
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